Evidence of local magnetic order in hep iron from Raman mode splitting 
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Experimental measurements of Raman spectra for hep iron at high pressure show two modes over 
a considerable pressure range in contrast to the prediction of one doubly degenerate mode for the 
hep lattice. We use density functional theory to investigate the influence of magnetic order on the 
Raman active modes of hep iron. We find an antiferromagnetic state that lifts the degeneracy of the 
transverse optical mode, and yields stable antiferromagnetic moments up to approximately 60 GPa 
(55 Bohr 3 ). The resulting frequencies of the two transverse optical modes are in good agreement 
with the experimental Raman shifts, lending support to the existence of local antiferromagnetic 
order in hep iron. 

PACS numbers: 75.30.-m,75.50.Bb,75.50.Ee,78.20.-m,78.20.Bh 



The presence of magnetism strongly influences mate- 
rial properties and phase relations in the 3d transition 
metals. This is particularily true for iron where the com- 
petition between magnetic and non-magnetic contribu- 
tions to the internal energy, differences in the vibrational 
and magnetic entropy, and differences in volumes all con- 
tribute to phase stability Q. Ferromagnetism stabilizes 
the body centered cubic structure over the close packed 
phases at ambient condition and an understanding 
of the spin-density wave ground state in its high tem- 
perature polymorph, cubic close packed (fee), has im- 
portant implications on its physical properties ||. Mag- 
netism in the high pressure polymorph, hexagonally close 
packed (hep) , is less well characterized and has long been 
thought to be absent, but its proposed presence [Q may 
play an important role in our understanding of material 
properties at pressure, with applications to the study of 
planetary interiors |Q and impact phenomena. 

The hep lattice has two atoms per unit cell giving rise 
to optical phonon modes through Brillouin zone folding. 
Advances in experimental techniques have made it possi- 
ble to measure the Raman active transverse optical (TO) 
mode in iron through optical spectroscopy || Q] . In con- 
trast to the fundamental prediction that this mode (E2 g ) 
be doubly degenerate, one set of experiments at high 
pressure show two peaks in the Raman spectrum up to 
pressures of 40 GPa 0. While amplitude and sharp- 
ness of the two peaks differ significantly, this observation 
nonetheless suggests that the symmetry of hep iron is 
lower than the atomic arrangement; spin-phonon inter- 
actions provide a symmetry-breaking mechanism. Based 
on first principles theory the lowest energy state found 
to date is an anti-ferromagnetic (afm) structure with or- 
thorhombic symmetry and four atoms in the unit cell 
(afmll) [Q. The afmll structure corresponds to a spin 
wave with wave vector q = (0, 1/2, 0), at the M point on 
the Brillouin zone boundary (Fig. [l]). This structure has 
been predicted to be stable for hep iron up to pressures of 
almost 60 GPa W . The magnetic structure results in two 



TO zone center modes (^2g) both of which are Raman 
active. 

Here we investigate the influence of afmll ordering 
on the zone center TO phonon mode over the com- 
pression range where finite moments arepredicted (>60 
Bohr 3 /atom). We calculate the zone center Raman fre- 
quencies by the frozen phonon approximation: energy 
changes are evaluated in response to small displacements 
along the phonon eigenvector, the second order term 
yielding the frequency. We base our analysis on to- 
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FIG. 1: Orthorhombic antiferromagnetic groundstate of hep 
iron (afmll). Open circles show the atomic positions at 2=1/4, 
filled circles at z~3/4 with the arrows indicating the direction 
of spin on the atoms. The orthorhombic unit cell is outlined 
with the axes given. The direction of the wavevector q for 
the afmll spinwave is along the 6-axis. The a- and &-axes also 
define the eigenvectors for the TO modes (TO a and TO{). 
The c-axis is out of the plane. 
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tal energies obtained with the full-potential linearized- 
augmented plane-wave method (LAPW) M using a gen- 
eralized gradient approximation |9| to the exchange corre- 
lation potential. We treat 3s, 3p, 3d, 4s, and 4p states as 
valence electrons for all volumes and use Rmt=2-0 Bohr 
for the muffin tin radii, RMrKmax — 9.0, a 12x6x12 
special k-point mesh, and a temperature broadening of 5 
mRy. With this set of computational parameters we have 
previously established convergence of relative energies to 
within 0.1 mRy/atom and magnetic moments to better 
than 0.05 /Us/atom [Q. The TO modes are characterized 
by displacements of the close packed planes with respect 
to one another. The eigenvectors for the the two A2 g 
modes in the afmll structure are along the orthorhombic 
a- and fr-axes (Fig. Q). For comparison we also compute 
the TO mode frequency for the non-magnetic case. To 
efficiently evaluate the energetics of the system we fixed 
the axial ratio in our calculations to c/a = 1.6, close to 
the experimentally (h], [tl] and theoretically Q estimated 
equilibrium value. 

Inspection of the afmll structure reveals that funda- 
mentally different spin interactions are involved in the 
two TO modes (Fig. |l|). For displacements along the or- 
thorhombic a-axis (TO a ) atoms approach nearest neigh- 
bors with unlike spin. For displacements along the or- 
thorhombic 6-axis (TOb) atoms alternately move towards 
a nearest neighbor with like spin and towards a pair of 
nearest neighbors of opposite spin. The resulting energy - 
displacement relations reflect the antiferromagnetic spin 
interactions (Fig. ^): for unlike spins approaching the 
energy is reduced with respect to the non-magnetic case, 
and for same spin nearest neighbor interactions (along 
positive b) there is an additional repulsion. The calcu- 
lated equilibrium magnetic moments are consistent with 
these findings (Fig. ^): They increase for displacements 
along a, as they do along negative b; for positive b they 
decrease considerably as the like spin nearest neighbors 
approach. 

Wc find that the TO a mode frequency agrees well with 
that of the lower frequency, higher amplitude peak found 
in the Raman experiments, and that the TOb mode fre- 
quency corresponds to the experimentally observed satel- 
lite peak at higher frequency (Fig. . The magnitude of 
the predicted TO mode splitting decreases as the afm 
moment is reduced by compression, in excellent agree- 
ment with the observation in the Raman experiments 
(Fig. ||). The magnitude of splitting is related to the di- 
rectionality of the spin wave considered here (with wave 
vector M), but other spin waves with wavevectors along 
high symmetry directions in the base of the hexagonal 
Brillouin zone will generally result in splitting of the TO 
mode. The systematic offset of the calculated frequencies 
by approximately 20 cm -1 (< 10%) is typical for a com- 
parison of computed and measured phonon frequencies 

Spin-phonon interactions have been found to have 



a strong effect on Raman scattering in a number of 
systems, including cupric oxide 1 14 and the copper- 
ruthenium oxide RuS^GdC^Os |15[ |. The general char- 
acter of the effect of spin-phonon interaction on the Ra- 
man spectra in these materials is consistent with the ob- 
servations for hep iron: broad, low amplitude, satellite 
peaks appear as the sample is cooled below the Curie or 
Neel temperature, upon further cooling such peaks gen- 
erally sharpen and gain amplitude as the degree of mag- 
netic ordering increases. Observation of an increase in 
amplitude of the satellite peak in hep iron upon cooling 
would lend support to the predicted magnetic structure. 

Other experimental investigations of possible magnetic 
states in hep iron have been inconclusive. Magnetism 
in hep irejii^x- ists at , over-qxpanded volumes in epitax i- 
ally grown multilayers on a ruthenium substrate Jl6[ al- 
though the na^re 5f-t^g spin ordering is still contr^ 

Nucle^X-r'ay^bsorption experirngffts |1| 
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FIG. 2: Energy of displacement of hep iron at V=70 
Bohr 3 /atom in the lower panel. Relative energy changes as 
a function of displacement for the non-magnetic structure 
(squares), and the afmll structure for displacements along 
TO a (triangles up) and TOb (triangles down). Dotted and 
solid lines are third order polynomial fits to the results. The 
upper panel shows the change in the afmll moment as a func- 
tion of displacement for TO a and TOb (same symbols, dashed- 
line to guide the eye). 
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FIG. 3: Raman frequencies as a function of atomic volume. 
Non-magnetic calculations are shown in filled squares. The 
afmll structure results in transverse optical frequencies (filled 
triangles) with TO a being the lower (up) and TOb (down) 
the upper branch. The solid lines through TO a and TOb 
are third order polynomial fits in V~ 2 ^ 3 . Experiments (Ref. 
^|) identify two peaks in the Raman spectra up to 40 GPa in 
open symbols. The stronger peak is shown in the triangles up, 
and the weaker peak in triangles down. The inset compares 
the split in Raman frequencies from theory (filled diamonds) 
and experiment (open diamond). 



be interpreted to show no moments in the high pressure 
polymorph: the change in absorption spectra is due to 
changes in the density of states as well as to spin related 
satellites. Diamond anvil cell in-situ Mossbauer measure- 
ments on hep iron have not shown any ordered magnetic 
ground state in its stability field |Eol 21, E2l E3I but it 



was recognized that magnetism cannot unambigously be 
ruled out based on the data ||(], The Mossbauer 

data suggest that correlation times of spin fluctuations 
must be relatively short, or that competing contributions 
to the effective field cancel each other as is the case in 
CaRuC>3 perovskite (24|, |25|. The in-situ observation of 
super-paramagnetism in hcp-Fe p2| may support the no- 
tion that the effective field is approximately cancelled, 
and suggests, following the procedure of Felner et al. p5| , 
a Mossbauer experiment on an appropriately doped hcp- 
Fe sample (e.g. with Ru). The inherent frustration of the 
triangular lattice might lead one to suspect more complex 
spin arrangments than those considered here, such as in- 
comensurate spin waves, as in the case of fee iron ^ p6| , 
non-colinear structures, spin glass, or a combination of 
these f27|. 

Acoustic modes can also be influenced by spin-phonon 
interactions as is evidenced by anomalous phonon disper- 



TABLE I: Raman shift for hep iron for the non-magnetic E^g 
mode and the two antiferromagnetic TO modes (A2 g : TO a , 
TOb) as a function of atomic volume. Pressure P is from the 
afmll equation of state (Ref. ^). 
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sion in fee iron near the zone center p8| . Indeed we find 
that the afmll ordering also influences the compressional 
properties of hep iron Q : the bulk modulus at pressure 
is reduced, bringing it in closer agreement with experi- 
mental data p0[ |j~H than non spin-polarized calculations. 
Better agreement of aggregate elastic properties between 
afmll and experimental data lends further, independent, 
support to the presence of magnetism in hep iron. 

The possible existence of magnetic correlations in hep 
iron may be important in our understanding of the recent 
observation of superconductivity in hep iron under pres- 
sure [^{J. Long thought to be antithetical, superconduc- 
tivity and magnetism have recently been observed simul- 
tanously in a number of systems |p0| , |3l| , |32| . Of particu- 
lar interest is the observation for ZrZn2 where magnetism 
and superconductivity appear to be directly coupled ]3lf| 
as is evidenced by the loss of superconductivity and mag- 
netism at the same pressure, leading to the speculation 
that hep iron might behave similarly J33|] . 

In conclusion, we report first-principles results on the 
zone center TO phonon frequency of hep iron for non- 
magnetic and afm structures at ambient and high pres- 
sure from a frozen phonon calculation. We find good 
agreement with experimental results for the absolute fre- 
quencies, and spin-phonon interactions provide a quan- 
titative explanation for the split of the TO mode in 
hep iron which is observed up to a pressure of 40 GPa. 
In combination with a considerably better agreement 
of elastic properties between experiment and the or- 
thorhombic afm structure, as compared to non-magnetic 
calculations, this suggests that magnetism may play an 
important role in the physical behavior of hep iron at 
high pressure. 
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